function [x_next, y_next, z_next] = SIMM(x, y, z, a, b, c)
    x_next = a * sin(b * z) * sin(c / x);
    y_next = a * sin(b * x_next) * sin(c / y);
    z_next = a * sin(b * y_next) * sin(c / z);
end